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Abstract 

We study the Collins effect in the azimuthal asymmetric distribution of hadrons inside a high 
energy jet in the single transverse polarized proton proton scattering. From the detailed analysis 
of one-gluon and two-gluon exchange diagrams contributions, the Collins function is found the 
same as that in the semi-inclusive deep inelastic scattering and e~^e~ annihilations. The eikonal 
propagators in these diagrams do not contribute to the phase needed for the Collins-type single 
spin asymmetry, and the universality is derived as a result of the Ward identity. We argue that this 
conclusion depends on the momentum flow of the exchanged gluon and the kinematic constraints 
in the fragmentation process, and is generic and model-independent. 
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I. INTRODUCTION 



Single-transverse spin asymmetries (SSA) in hadronic processes have a long history 
Recent experimental measurements of SSAs in polarized semi-inclusive lepton-nucleon deep 
inelastic scattering (SIDIS) [3|, in hadronic collisions , and in the relevant e~^e 

;l 

annihilation process [Si], have renewed the theoretical interest in SSAs and in understanding 
their roles in hadron structure and Quantum Chromodynamics (QCD). There are several 
approaches to understanding SSAs within the QCD framework |9|, [lO|, lul- Transverse- 
momentum-dependent (TMD) parton distributions and fragmentation functions, and their 
relevance for semi-inclusive DIS, the Drell-Yan process, di-hadron production in e^e~ annihi- 



ations, anc 
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the single-inclusive hadronproduction at hadron colliders have been investigated 
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19|, 
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26]. 



Two important contributions from these TMD parton distributions and fragmentation 
functions have been mostly discussed in the last few years: the Sivers quark distribution 



and the Collins fragmentation function. The Sivers quark distribution 16|] represents a dis- 
tribution of unpolarized quarks in a transversely polarized nucleon, through a correlation 
between the quark's transverse momentum and the nucleon polarization vector. The exis- 
tence of the Sivers function requires final/initial-state interactions 20|], and an interference 



between different helicity Fock states of the nucleon [20|, 127[. The Collins function repre- 
sents a correlation between the transverse spin of the fragmenting quark and the transverse 
momentum of the hadron relative to the "jet axis" in the fragmentation process. Like the 
Sivers function, it vanishes when integrated over all transverse momentum. 

One of the most nontrivial properties associated with the Sivers and Collins functions 
are their universality properties. Although they both belong to the so-called "naive-time- 
reversal-odd" functions, they do have different universality properties. For the quark Sivers 
function, because of the initial/final state interaction difference, they differ by signs for the 
SIDIS and Drell-Yan processes 2^, [21, 28|. This non-universality has also been extended to 
other processes, such as the dijet-correlation in hadronic reactions, where it was found that 
both initial and final state interactions contribute to the SSA, and there exists non-trivial 
relation between this and those in the SIDIS and Drell-Yan processes 



standard TMD factorization breaks down 
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On the other hand, there have been several studies showing that the Collins function 
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is universal between different processes, primarily in the SIDIS and e^e annihilation 32, 



33 



34| . In these discussions, the gauge links in the fragmentation functions do not play a 
crucial role to leading to a nonzero Collins function though they are important to retain the 
gauge invariance, whereas it has been well understood that the gauge links in the parton 
distributions play very important roles to obtain non-zero quark Sivers function. 

The Collins effect in the fragmentation process and its universality has been recently 



extended to the hadron production in pp collisions 



351], where the azimuthal distribution of 



hadrons inside a high energy jet can probe the Collins fragmentation function and the quark 



transversity distribution [36|] in the single transverse polarized nucleon-nucleon scattering. 
In this paper, we will give the detailed derivation of these results, and argue that the 
universality is in general and model- independent. 

We are interested in the hadron production from the fragmentation of a transversely 
polarized quark which inherit transverse spin from the incident nucleon through transverse 



spin transfer in the hard partonic scattering processes 37|, l38|, l39|]. As shown in Fig. 1, we 
study the process, 

p{Pa, S±) + p{Pb) ^ jet{Pj) H{Ph) + X , (1) 

where a transversely polarized proton with momentum P4 scatters on another proton with 
momentum Pb, and produces a jet with momentum Pj (transverse momentum P± and 
rapidity yi in the Lab frame). The three momenta of Pa-, Pb and Pj form the so-called 
reaction plane. Inside the produced jet, the hadrons are distributed around the jet axes. 
A particular hadron H will carry certain longitudinal momentum fraction of the jet, 
and its transverse momentum PhT relative to the jet axis will define an azimuthal angle 
with the reaction plane: (phi shown in Fig. 1. Thus, the hadron's momentum is defined 
as Ph = ZhPj + PhT- The relative transverse momentum PhT is orthogonal to the jet's 
momentum Pj: PhT ■ Pj = 0- Similarly, we can define the azimuthal angle of the transverse 
polarization vector of the incident polarized proton: (pg. 

The leading order contribution to the jet production in pp collision comes from 2-^2 
sub-processes, where two jets are produced back-to-back in the transverse plane. For the 
reaction process of (1), one of the two jets shall fragment into the final observed hadron. 
In this paper, we study the physics in the kinematic region of PhT ^ -Pi- The unpolarized 
cross section contribution from the partonic 2^2 process ab —>■ qc where the final state 
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FIG. 1: Illustration of the kinematics for the azimuthal distribution of hadrons inside a jet in pp 
scattering. 



quark q fragments into final observed hadron H, can be written as 

b=q,9 

where dV.S. = dyidy2dP]_dzd'^ P^t represents the phase space for this process, yi and y2 
are rapidities for the jet Pj and the balancing jet, respectively, P± is the jet transverse 
momentum, and the final observed hadron's kinematic variables Zh and PhT are defined 
above. Here, x and x' are the momentum fractions carried by the parton "a" and "6" from 
the incident hadrons, respectively. In the above equation, fa and fb are the associated 
parton distributions, and Dq{zh, Pht) is the TMD quark fragmentation function. The hard 
factors Hab^qc are equal to the partonic differential cross section for the relevant subpro- 
cess: Hab-*qc = da / di\ab^qc- Similarly, the differential cross section for the transverse-spin 
dependent scattering process can be written as 

da{S^) _ J2^'f^^^r^^SqT{x)5q{z,,P,Tf^'- 



dV.S. ^ ' ' ' ^ ' ' ' ' Mh 

b=q,g 



X 



p/3 _ Pb ■ PhT pl3 



X J , (3) 

= e'^'^'^PAt.PBu/PA ■ Pb with convention e^^^s = and H^b^^b is the hard factor for the 
partonic channel qb — * qb. Here, 6qT{x) (also noted by 6q, hig and A^q in the literature) 
is the quark transversity distribution, and 6q the Collins fragmentation function 17|] (also 
noted as AD or in the literature). 

It was argued that the Collins function is universal between the above process and other 
processes such as e~^e~ annihilation and SIDIS [sS^. As an example, we will demonstrate 
this universality for the particular partonic channel qq' qq' contribution to our process, 
and all other channels will follow accordingly. For convenience, we list the hard factors for 
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(a) 



(6) 



^^'IG. 2: Quark fragmentation to pion production (a), and in pp scattering in a model described in 
this channel, 

for the unpolarized and single-transverse-spin polarized cross sections, respectively. Here s, 
t, and u are the usual partonic Mandelstam variables. 

The rest of this paper is organized as follows. In Sec. II, we calculate the differential cross 
sections for the unpolarized and single-spin dependent scattering processes from qq' qq' 
channel contributions, and demonstrate the universality of the Collins function. Especially, 
we will present a detailed calculation for one-gluon exchange diagrams which are essential 
for the universality argument. An extension to two-gluon exchange diagrams is presented in 
Sec. III. We summarize our paper in Sec. IV. 



II. UNIVERSALITY OF THE COLLINS FUNCTION 



We follow the model used in Ref. 17|] to calculate the quark fragmentation into a pion. As 
shown in Fig. 2(a), a quark (with momentum k) fragments into a pion (with momentum Ph) 
by the vertex from a model described in 40| . A simple calculation will give the unpolarized 
quark fragmentation function 17l |. 



(5) 



where g is the coupling between the quark and pion, M is the quark mass. 

We can also use this model to calculate pion production in hadronic process of (1). In 
Fig. 2(b), we show the Feynman diagram for the typical partonic channel qq' — ^ qq' contribu- 
tion, where the initial quarks have momenta Pa and Pb, respectively. In the final state, the 



produced pion has momentum Ph, and the associated final state quark has momentum k', 
whereas the balancing jet has momentum P2. We further introduce a light- like momentum 
k: k^ = 0, which represents the dominant component of the fragmenting quark's momentum. 
It can be parameterized as follows, 

k^~PA-iPB + P±, (6) 
s s 

where P± is the transverse momentum for the fragmenting quark in the Lab frame, s, t and 
u as mentioned above, are the usual partonic Madelstam variables for this partonic process: 
s — 2Pa ■ Pb: i — —2Pa ■ k, and u — —2Pb ■ k. In our discussions, the jet's transverse 
momentum (in the lab frame) is the large momentum scale at the same order as s, t 
and u. Of course, the full momentum of the fragmenting quark Pi = Ph + k' is off-shell in 
this diagram. However, its off-shellness is much smaller than Pj_. In order to formulate the 
final state hadron's momentum, we introduce a conjugate light-like vector h: n° = k^ and 
n — —k. It is convenient to define this momentum in the center of mass frame of the two 
incident momenta Pa and Pb- In this frame, we have 

n^P,^-ipA-^PB-P± , (7) 
s s 

which happens to be the momentum of the balancing jet. From above, we have A;^ = = 
and k ■ n — s/2. In the following calculations, we will work in this particular frame. We 
emphasize that our results do not depend on the frame. 

With the above two momenta, we can formula the final state pion's momentum as 

Ph = Zhk + + PhT , (8) 

2zhk ■ n 

where Zh — Ph ■ n/k ■ n is the momentum fraction of the fragmenting quark carried by the 
pion in the final state, PhT is the transverse momentum relative to the fragmenting quark 
momentum k: PhT • k — and PhT ■ fi — 0. In the above parameterization, we have neglect 
the pion mass, which is not relevant in our calculations. Similarly, we can formulate the 
associated final state quark momentum k' as, 

k'^{l- Zh)k + -P^-^n - PhT , (9) 
2(1 — Zh)k ■ n 

where we have kept the quark mass, because it will be relevant for the nonzero single spin 
asymmetry discussed below. 
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(a) 



I k. sx 




■H^ /"'\ 



(d) 



FIG. 3: Universality of the Collins function in e~^e~ (a), deep inelastic scattering (h), and pp 
scattering (c), when we have dressed quark propagator associated with the fragmenting quark in 
these processes. The universal Collins function can he calculated from the diagram in (d). The 
blobs in the diagrams represent the dressed quark propagator in this model. 



With the above decompositions for the relevant momenta, it is straightforward to calcu- 
late the Feynman diagrams for this process in Fig. 2(b). In the calculations, we will utilize 
the power counting method to keep the leading order contributions, and neglect all higher 
order corrections of Pht/Pl or MjP^. By doing that, we can separate the short distance 
physics (at the scale of Pi) from the long distance physics (at the scale of Put and M). 

Finally, the cross section contribution from Fig. 2(b) will be. 



da'' 



12fs2 



u 



dV.S. 



Am 



167r2 + zlM^ ' 



(10) 



in the limit of Phr <^ Pi. This result is indeed factorized into the hard factor H':^3_„„, in 



qq'^qq' 



Eq. (jll) times the fragmentation function in Eq. ([5]) calculated from Fig. 2(a). 



Now, we turn to discuss the SSA in this process. We need to generate a phase from the 
scattering amplitudes to have a non- vanishing SSA. As suggested in [17], the dressed quark 
propagator in this model may contribute to such a phase. Similarly, the vertex correction to 
the quark-pion vertex can also contribute a phase 411] • If the phase comes from the above 
sources, it is easy to argue the universality of the Collins function between our process and 
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the SIDIS/e"'"e~ process, because they are the same. For example, as we show in Fig. 3, the 
dressed quark propagator associated with the fragmenting quark can contribute to a nonzero 
phase [l7| , which will contribute the same to the Collins function in all these three processes: 
e^e~ annihilation, SIDIS, and hadron production in pp scattering. This propagator can be 



parameterized as: i{A /Pi + BM)/{Pf — M^) 17|, where A and B are complex numbers. 
Following the above calculations for the unpolarized cross section, we will find the single-spin 
dependent cross section for process (1) from Fig. 3(c) can be written as 



p/3 _ Pb ■ PhT pf3 

Pb-Pj ' 



oI-kNI-IAsu zljl- Zh)2Mlm{A*B) 

where again we only keep the leading order contribution in the limit of Pht ^ -Pi and 
M <C Pj_. In the derivation of the above result, the following identity has been used to 
simplify the final expression. 



sP± ■ S^P^P^T +{u- t)PB ■ PhxPlSi - iuP^rSi = , (12) 



which holds in our working frame. The above differential cross section can be factorized into 
the Collins function calculated from the dressed quark propagator from Fig. 3(d) [l?! and 
the hard factor from Eq. (jlj) in this partonic channel qq' —>■ qq', and this Collins 

function will be the same as that in e~^e~ and SIDIS processes in Fig. 3(a) and 3(b). 

Similarly, the vertex corrections contributions to the Collins function can be analyzed 
accordingly, and the same factorization and universality of the Collins function will follow. 

The main issue of the universality discussion concerns the extra gluon exchange contri- 
bution between the spectator and hard partonic part jQ. For example, in our case, because 
the hadron is colorless while the quark is colored, the remanet in the fragmentation process 
will be also colored. Thus the gluon exchanges between the remanet and the other parts of 
the scattering amplitudes become essential. In Fig. 4, we have shown all these interactions, 
including the gluon attachments to the incident quarks (a,c) , and final state balancing quark 



(d) and the internal 
those discussed in 



uon propagator (b). These diagrams are much more complicated than 
32| for SIDIS and e~^e~ processes, where there is only one diagram con- 
tribution in both cases. Therefore, the universality argument for the Collins function is not 
straightforward. However, the dominant contribution to the fragmentation function comes 
rom the kinematic region where the exchanged gluon is parallel to the final state hadron 



42]. Otherwise, their contributions will be power suppressed in the limit of PhT ^ -Pi or 





FIG. 4: Gluon exchange diagrams contributions to the Collins asymmetry in pp collisions. The 
short bars indicate the pole contributions to the phase needed for a non-vanishing SSA. The addi- 
tional two cuts in (d) cancel out each other. 

belong to a soft factor. For these collinear gluon interactions, we can use eikonal approxi- 
mation and Ward identity to sum them together to form the gauge link in the definition of 
the fragmentation function |42 |. 



Meanwhile, we notice that the contributing phases of the diagrams in Fig. 4 come 



rom 



the cuts through the internal propagators in the partonic scattering amplitudes [20|, . In 
Fig. 4, we labeled these cut-poles by short bars in the diagrams. From our calculations, we 
find that all these poles come from a cut through the exchanged gluon and the fragmenting 
quark in each diagram, and all other contributions either vanish in the leading order con- 
tribution or cancel out each other. For example, in Fig. 2(d), we show two additional cuts, 
which contribute however opposite to each other and cancel out completely. To see this 
cancellation more clearly, we can write down the momentum integral of the exchange-gluon. 



(4^ 

X 



M{q) 



1 



1 



1 



{k' - g)2 _ M2 + ie g2 + ie {P^ + qf + ie 
1 1 



(P2 -Pb + qf + le (Ph + k' - qf - + le ' ^^^^ 
where M.{q) represents the denominators coming from the scattering amplitude. By power 
counting analysis, the dominant contribution to the fragmentation comes from the kinematic 



region of q being parallel to the final state hadron's momentum q ~ Ph. Prom this fact, we 
can parameterize q in terms of k and n, and define g+ — q • fi/k • fi and q" — q ■ k/k • fi. 
Thus the integral of momentum q becomes d'^q — k-ndq~^dq~d'^qT, where qr is the transverse 
momentum relative to the jet momentum k. Because q is parallel to k, q^ will be order 1, 
whereas q~ will be order of q^/Q^- When wc perform the integrals of dq^dq~ , we need to 
take two poles from the above propagators to obtain a nonzero Collins asymmetry. These 
poles will form a cut through the Fenyman diagram. Physically, these cuts represent the 
kinematic allowed final state re-scattering in the diagram. 

By examining the behaviors of the propagators in the above kinematic region, we further 
notice that the t— channel gluon propagator 1 / (P2 ~ Pb + does not contribute to a pole. 
This is because this propagator is far off-shell: (P2 — PbY = t ~ — I-PlP- If we take a 
pole from this propagator, we have to constrain the momentum of q being proportional 
to P2 and Pb-i whose contribution will be power suppressed. Thus, we shall calculate the 
pole contributions from other propagators. In Fig. 4d, we show three possible cuts which 
are kinematic allowed for this diagram. Two of them are associated with the propagator 
l/(-P2 + lY- This propagator involves large momentum P2, and can be simplified by using 
the eikonal approximation, 

111 
{P2 + qy + ie^2P2-kq^ + ie' 

The pole contribution from this propagator is proportional to S{q'^). With this delta func- 
tion, the integral over q~ vanishes, because the rest poles are in the same half plane of 

dq- 1 1 



I 



271 {k' - qf -M^ + ie {Ph + k' - qf - + ie 

I I = (15) 

J 27r 2{k'+-q+)q- + ■■■ + ie2{P+ + k'+- q+)q- + -- - + 16 ' ^ ^ 

where we have used the fact that — g+ > and -|- — g"^ > 0. This means that 
the two cuts associated with the propagator 1/(P2 + q)"^ cancel out each other. The above 
result depends on the momentum flow of q in this diagram and the time-like process in the 
fragmentation region requiring that k'^ > and P^ > 0. 

Therefore, the only contribution to the nonzero SSA associated with the Collins effect 
comes from the cut going through the fragmenting quark and the exchange-gluon, as we 
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labeled by short bars in this diagram. Summarizing the above analysis, we find that the 
contribution from this diagram can be written as 

k-h r dq^dq~cPqT 11 1 

(P2 - hY{P2 - PbY J ^l-g+ (A;'-g)2-M2-^^^^ 

x5{q^)5{{P^ + k'-qf-M^) , (16) 

where we have made the eikonal approximation for the propagators 1 / (P2 + 0')^ and 1 / (P2 ~ 

Similar analysis can be done for all other diagrams, and we find that their contributions 
come from the same poles of the fragmenting quark and the exchange-gluon. Therefore, their 
contributions will have the similar expression as Eq. ( fT6l) with the same delta functions in 
the integral: S{q'^)6 {{Ph + k' — g)^ — M^) and the propagator l/{{k' — g)^ — M^). Thus the 
contributions from all these diagrams can be summed together. In this sum, we notice that 
the different diagrams have different color-factors, 

• ^Tr [T^T^T^T'^] x Tr[T"T^] = Cp x + -L^^^Tr [T^T'T^] , 

4(b) : -^TriT'^T^T^] x Tr[T'^T'^]^/,,, = _i^!:^Tr[T"T^T'=] , 

4(c) : -^Tr [T"T^T*] x Tr [T'^T^T"] , 

c 

4(d) : ArTrfT'^T'^T^] x Tr[T'^T"T^] . (17) 

We further find that the contributions (without the color-factors) from the diagrams (c) and 
(d) are opposite to each other. Thus, their total contribution will be the difference on the 
color-factor, which is -^^'^Tr[T"T^T'^]. That means the contributions from all these four 
diagrams can be grouped into two terms with different color factors: one with Cp x 
and one with ^%^^Tr[T"T''T'=]. The latter one vanishes in the leading order of PhT/P± 
after we sum all diagrams contributions, and thus we are left with the first color-factor 
contribution. 

After summing over all diagrams' contribution, the spin-dependent differential cross sec- 
tion coming from the Collins effect will be 

da{S^) Nl - 1 Uu 



dV.S. TT Am -p 



Pb-Pj ' 



g"^ ^ 2 [ dq^dq d'^qr 
(277)3 P9sJ (27r)4 



{k'-q)^ 
11 



where Qs is the strong couphng. From the above result, we find a clear separation of the 
short distance physics at the scale Pi and long distance physics at the scale PhT- The 
short distance part is just the hard factor H^^^}^^^, for the spin- dependent cross section, 
which can be calculated from the partonic process with both initial and final state quarks 
transversely polarized 38|, as we show in the left panel of Fig. 5. The long distance part of 
the above result can be factorized into the Collins fragmentation function calculated from 
the right panel of Fig. 5. In this part, because the integral is proportional to P^j,, we can 
combine the two terms in the integral into one expression contained in the Collins function. 
Therefore, the spin-dependent cross section Eq. ( |T8l) can be re-written as 



dV.S. vr Am -£2 ^ 



,/3 Pb ■ PhT „/3 



where the Collins function 5q is calculated from the Feynman diagram in the right panel of 
Fig. 5, 

, p ^ g"" 2^ [ dq+dq-(fqT ( + g „^ 

ZH,Pht) - pa^J^39sCFj (2vr)4 PkT-ZkqT) 

^ -S{q')6{{P, + k'-q)'-M') , (20) 



{k' - qf - M2 

where the index a is not understood as a sum. This final result demonstrates that we do 
have a factorization for the spin-dependent cross section into the hard factor H^^^}^^^, times 
the Collins fragmentation function, and the Collins function is the same as that in e^e~ and 
SIDIS processes js^. 

Therefore, by using the Ward identity at this particular order, the final results for all 
the diagrams of Fig. 4 will sum up together into a factorized form as shown in Fig. 5, 
where the cross section is written as the hard partonic cross section for q{Sj_)q' q{sj_)q' 
subprocess multiplied by a Collins fragmentation function. The exchanged gluon in Fig. 4 



14l | as shown in 



is now attaching to a gauge link from the fragmentation function definition 
the right panel of Fig. 5. 

The key steps in the above derivation are the eikonal approximation and the Ward iden- 
tity. The eikonal approximation is valid when we calculate the leading power contributions 
in the limit of PhT <^ kj_. The Ward identity ensure that when we sum up the diagrams with 
all possible gluon attachments we shall get the eikonal propagator from the gauge link in 
the definition of the fragmentation function. The most important point to apply the Ward 
identity in the above analysis is that the eikonal propagator does not contribute to the phase 
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FIG. 5: Factorize the contributions from Fig. into the hard partonic cross section multiplied by 
the universal Collins fragmentation function. The short bars indicate the pole contribution to the 
Collins function. 

needed to generate a nonzero SSA. This is what we have shown for the Colhns asymmetry in 
the above calculations, and the reason, as we mentioned above, is due to the momentum flow 
of the exchanged gluon and the kinematic constraints in the fragmentation process. We will 
show in the next section, that for the two-gluon exchange diagrams the eikonal propagators 
do not contribute to the phase for the nonzero Collins SSA in this process. Therefore, we 
conjecture that the above conclusions are valid to higher order contributions too. 

This argument can not apply to the SSA associated with the parton distributions, where 
the eikonal propagator does contribute to the phase to generate a nonzero SSA. That is the 
reason we have sign differences for the Sivers functions in SIDIS and Drell-Yan processes. 



III. TWO-GLUON EXCHANGE CONTRIBUTIONS 

As we discussed in the last section, to demonstrate the universality of the Collins function, 
we have to apply the Ward identity to sum up all gluon exchange contributions into the 
gauge link from the definition of the fragmentation function. In order to use this argument, 
the eikonal propagator should not contribute to the phase needed to generate nonzero SSA 
associated with the Collins effects. This has been explicitly demonstrated in the last section 
for the one-gluon exchange contribution. In this section, we will extend the discussions to the 
two-gluon exchange contributions. Especially, we will show that these eikonal propagators 
do not contribute to the phase for the SSAs. The reason, again, is due to the time-like 
feature and the momentum flow in the fragmentation process. 

We will focus our discussions on some representative diagrams from the two-gluon ex- 
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FIG. 6: Example diagrams for two-gluon exchange contributions (a,h,c); and one real gluon radia- 
tion contributions (d,e,f). 

change contributions. All other diagrams will follow accordingly. We show these diagrams 
in Figs. 6(a,b,c). The contribution from Fig. 6(a) will depend on the following integral of 
the exchange gluons' momenta qi and g2, 

fd^d^j^. . 1 1 1 

J (27r)4 (27r)4 (P^ - q,y + le {Pa - qi - gs)^ + ^e {k' - q^f + le 

1 111 
(21) 

{k' - gi - q2y + ie{k - qi- ^2)^ + gf + ie ql + ie 

where k = Pi = k' + is the fragmenting quark's momentum and A^(gi,g2) represents 
the numerators depending qi and 52? especially their transverse momentum components. 
Following the arguments used in the last section, the first two propagators in the above 
expression can be further simplified by using the eikonal approximation, and then we will 
obtain the following expression 

dq^dq^ dq^dq^ 111 1 
(27r)2 (27r)2 -q+ ie -q+ - q+ + it {k' - qiY + ie {k' - qi - q2Y + ie 
111 

(22) 



{k - qi - q2Y + ieql + ie ql + ie 
where qf follow the definitions in the last section. The normalization of the above integral 
has been changed for convenience. This normalization is not relevant for our discussions, 
because we want to show that the eikonal propagators do not contribute to the phase needed 
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for a nonzero SSA, not the actual contribution from this diagram. We will show if we take 
pole contributions from these two eikonal propagators, the final integral will vanish. Because 
of the existence of two eikonal propagators, the analysis will be more complicated than that 
in the last section. We discuss their contributions separately. 

1. pole contribution from 



-<?! -<l2 

If we take pole of this eikonal propagator, and will be constrained: Qi +q2 = 0, 
and the integral of (l22l) will become, 

dqfdqt 5{qt + qt) f dq^ dq^ 1 1 



27r g+ J (27r)2 _2(A;'+ - + Ai + ze -2A;'+(gf + ^2") + A2 + ze 

1 11 

(23) 



-2k+{q^ + 52 ) + A3 + ze 2q{q^ + A4 + 2g^g^ + A5 + ie 

where Aj are some quantities depending on the transverse momenta of g^, k' and Ph- 
The following analysis does not depend on the details of these numbers. In deriving the 
above equation, we have used the constraint of + = to simplify the expression. 
We further notice that k'^ > and A;+ > 0. Thus, the poles of the second and third 
factor in the integral of gjf and are both in the upper half plane. If g^*" > 0, which 
means that g^ < 0, the pole of the fifth factor will be also in the upper half plane of 
g^. Therefore, the poles of the three factors (the second, third and fifth) depending 
on g^ are all in the upper half plane of g^, and the integral over g^ will vanish, and 
so will the above integral. Similarly, if qf < 0, the pole of the fourth factor will be in 
the upper half plane of gjf . Meanwhile, we will also have k'^ — q^ > 0, and the pole 
of the first factor will be in the upper half plane too. Therefore, the poles of the four 
factors (the first, second, third and fourth) depending on gf are all in the upper half 
plane of gf. The integral over g]" will vanish, and so will the above expression. In 
conclusion, in any case of gj,^ > or g^^ < 0, the above integral vanishes, and we do 



not have contribution from the pole of — + ^ . — -. 
2. pole contribution from 



-<?l +«<: 

Because gj*" = 0, we can simplify the integral of (!22l) as follows, 

dqfdq2 ^{(li) f dq^dq2 1 1 



271 g+ J (27r)2 -2A;'+gf + Ai + ze -2(A;'+ - g+)(gf + g2") + A2 + ze 
1 1 



-2{k+ - q+){q^ + q^) + A3 + ze 2g2+g2" + A5 + 
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(24) 



Again, the normahzation has been changed for convenience. Because > k'~^ , we 
will analyze the contributions of the above equation by classifying the different regions 
of (a) q2 > k"^', (b) < k'^; (c) k'^ < q^ < k^ . In the region of (a), we will 
have fc"^ — < and k'^ ~ it < 0- Therefore, the poles of the three factors (the 
second, third and fourth) are all in the lower half plane of q2 ■, and the integral over 
vanishes. In the region of (b), we have k'^ > 0, k'~^ — > and k^ — q^ > 0. Thus, 
the poles of the three factors (the first, second and third) depending on q^ are all in 
the upper half plane, and the integral over vanishes. In the region of (c), we have 
q2 > 0, k'^ — g^ < and k'^ — q^ > 0. Therefore, the g^ integral will pick up the pole 
of the third factor, which actually determines the value of gf + g2^. After substituting 
this back into the equation, we will find the second factor does not depend on q^ any 
more. The only dependence comes from the first factor. Obviously, this integral over 
q^ will vanish. In conclusion, in any case of (a,b,c), the above integral vanishes, and 
there is no contribution from the pole of — . 

-q{+ie 

In summary, there is no contribution to the SSA from the pole of the eikonal propagators 

in the diagram of Fig. 6(a). Similarly, the contribution from Fig. 6(b) depends on the 

following integral, 

dq^dqi dq^dq^ 111 1 
(27r)2 (27r)2 -q+ ie -q+ - q+ + ie {k' - q2Y + it {k' - qi - q2Y + it 

' ' (25) 



[k - qi - q2Y + itql + it ql + it 
where we have made the eikonal approximations for the two propagators along the incident 
quark line Pa- Comparing with Eq. (1221) . we find the only difference is the third factor 
Qi ~^ (l2- Again, we can show that none of the two eikonal propagators will contribute to 
the phase needed for a nonzero SSA. We will discuss their contributions separately. 

1. pole contribution from — + . — -. 

After taking this pole, the integral of ( l25l) will become, 

dqtdqt 5{qt + qt) f dq^dq2 1 1 



27r g+ J (27r)2 -2{k'+ - q+)q2 + Ai + it -2k'+{q^ + q2) + A2 + it 

(26) 



-2k+{q^ + g2 ) + A3 + ze 2gf g^ + A4 + ie 2g^g2 + A5 + 
which is the same as Eq. (1231) if we interchange gf and qf. Thus, the above integral 
will vanish by the same arguments we have used for Eq. (|23|1 . 
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2. pole contribution from — r — - 

Because = 0, we can simplify the integral of fl25|) as follows, 

dqfdq2 S{qi) f dqidq2 1 1 



2n q+ J (27r)2 -2{k'+ - q^)q2 + + ie -2{k'+ - q+){q^ + q^) + A2 + ie 

(27) 



-2{k+ - q+){q^ + q^) + A3 + ie 2^2+ ^2" + A5 + ' 
Again, we classify three different regions of q^ in the above equation: (a) q2 > k'^', (b) 
q2 < k'~^; (c) k'^ < < A;+. The contributions from (a) and (b) regions vanish by 
the same reasons as we have shown for Eq. in the above. In region (c), we have 

~ < and k~^ — q2 > 0. Therefore, the gf integral will pick up the pole of the 
third factor, which again actually determines the value of gjf + g^. After substituting 
this back into the equation, we will find the second factor does not depend on g^ any 
more. The only dependence comes from the first and last factors. Obviously, this 
integral over g^ vanishes because k''^ — g^ < and g^ > 0, and the poles of these 
two factors are both in the lower half plane. In conclusion, in any case of (a,b,c), the 
above integral vanishes, and there is no contribution from the pole of — 4 — . 

Similarly, the contribution from Fig. 6(c) will depend on the following integral, 

dq^dq^ dq^dq^ 111 1 
(27r)2 (27r)2 -q+ + q+ + [k' - qif + ie {k' - qi - q2y + it 
111 
(A; - gi - q2Y + le qf + le g| + le 

where again we have made the eikonal approximations. There are two eikonal propagators 
in the above equation, and as above we will discuss their contributions separately. 

1. the pole contribution from — t — 

After taking this pole, gj*" = 0, the above equation Eq. (125]) will reduce to Eq. ^IM . 
According to the same arguments we used there, there will be no contributions from 
this pole. 

2. the pole contribution from — 

This pole contribution means that g^ = 0, and the integral of Eq. (!28l) become 

dqt dqt 5{q^) f dq];dq2 1 1 

27r g+ J (27r)2 _2(A;'+ - g+)gf + Ai + -2(fc'+ - g+)(gr + gs") + A2 + 
1 1 



-2{k+ - qt){qi + g2 ) + A3 + ie 2q+q^ + A5 + ze 
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(29) 



which will be identical to Eq. (127|) if we interchange qf to q^. Using the same argu- 
ments there, the above integral vanishes. 

The above three examples are typical diagrams we encounter for the two-gluon exchange 
contributions for this channel. All these diagrams can be analyzed by a similar manner, 
and we will find that the eikonal propagators do not contribute to the phase needed to 
a nonzero SSA. Because of this fact, all these diagrams can be summed together to form 
the contributions from the gauge link in the fragmentation function, where the two gluons 
attach to the gauge link similar to the diagram we have shown in Fig. 5. Since there 
is no contributions from these eikonal propagators, the Collins function calculated from 
these diagrams will be the same as that in e~^e~ and SIDIS processes, and the universality 
preserved. 

We have also drawn some other diagrams at this order in Fig. 6(d,e,f), which contribute 
to a real gluon radiation in addition to the gluon exchange. The analysis of these diagrams 
also show that we do not get contribution from the pole of the eikonal propagators. For 
example, the contribution from Fig. 6(d) depends on 

rf^g 1 1 1 11 



(27r)4 -q+ + ie -q+ - kf + ie {k' -qf + it{k-q- kif + ie q^ + ie ' ^^^"^ 
where ki is the momentum for the radiated gluon. We have two eikonal propagators in the 
above equation. However, none of them contributes to the phase needed to a nonzero SSA. 

1. the pole contribution from — 

This pole contribution means that = 0, and the integral of g~ will reduce to 

/dq~ 1 1 

2tt -2k'+q- + Ai + ie -2{k+ - kt)q- + A2 + ie ' ^ ^ 

Because k^ = k^ + k'^ + > kf and k'^ > 0, the poles of the above two factors are 
both in the lower half plane of g~, and the integral vanishes. 

2. the pole contribution from — . \. . 

After taking this pole, we will have the following q^ integral 

J 2tx -2{k'+ - q+)q- + Ai + ie -2k+q- + A2 + ie 2q+q- + A3 + ie ' ^ ^ 
Because the pole constrains that q^ = —kf < and k'~^ — g+ > 0, the poles of the 
above three factors are all in the lower half plane. The integral over q~ vanishes. 
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In summary, there is no contribution from the pole of the eikonal propagators in the diagram 
of Fig. 6(d). 

The contribution from Fig. 6(e) will depend on the following integral, 

f d'<l ^ 1 I 1 1 .33) 

J (27r)4 -kf + ie -q+ - kf + ie {k' -qf + ie{k-q- k^f + ieq^ + ie ' ^ ' 

Because k'l > 0, we only have one possible pole contribution from the eikonal propagator 
l/{—q'^ — kf + ie), which vanishes by the same reason as above for diagram Fig. 6(d). 
Similarly, if the gluon with momentum q attaches to the radiated gluon instead of the 
incident quark line with momentum P4 (we did not show this diagram in Fig. 6), the 
contribution vanishes by the same reason. 

The contribution from Fig. 6(f) depends on the following integral 

d^q 1 1 1 11 



/ 



(34) 



(27r)4 -q+ + ie {k' - q^ + ie {k - q - kif + ie{k- qf + ieq^ + ie ' 

after eikonal approximation. If we take the pole contribution from the eikonal propagator, 
the above integral will reduce to 

dq- 1 1 1 



/ 



(35) 



27r -2k'+q- + Ai + ie -2{k+ - k^)q- + A2 + ie -2k+q- + A3 + ie ' 

Again, because k'^ > kf, the poles of the above three factors are all in the lower half plane, 
and the integral over q^ vanishes. Thus, there is no contribution from the pole of the eikonal 
propagator for this diagram. 

In summary, for the gluon radiation diagrams, there is no contributions from the pole of 
the eikonal propagators. Because of this fact, we can use Ward identity to sum all these 
diagrams together to form the gauge link contribution from the fragmentation function, 
similar to the diagram in Fig. 5 with an additional gluon radiation. 

Concluding the analysis of the two-gluon exchange diagrams in Fig. 6, the eikonal prop- 
agators do not contribute to the phase needed for the nonzero SSA associated with the 
Collins effect. Therefore, we can apply the Ward identity at this order to sum all these 
diagrams plus other similar ones. This sum will lead to the gauge link contribution from the 
fragmentation function definition, and the fragmentation function will be the same as that 
in e''"e~ and SIDIS processes. 
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IV. SUMMARY AND DISCUSSIONS 



In this paper, we have shown that the Colhns function in hadron production in single- 
transverse-spin polarized pp scattering is the same as that in e+e~ and SIDIS processes. 
This universality is a general and model-independent observation, and depends on the fact 
that the eikonal propagators do not contribute to the phase needed for a nonzero SSA. 
We have demonstrated this by explicit calculations for one-gluon exchange diagrams which 
corresponds to one eikonal propagator in the amplitudes, and two-gluon exchange diagrams 
which correspond to two eikonal propagators. Although our calculations were based on a 



model 
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40[ |. the analysis and arguments are quite general. The results, as we emphasized. 



depend on the momentum flow and kinematic constraints in the fragmentation process. 

This observation is very different from the SSAs associated with the parton distributions, 
where the eikonal propagators from the gauge link in the parton distribution definition play 
very important role. It is the pole of these eikonal propagators contribute to the phase 
needed for a nonzero SSA associated with the naive-time-reversal-odd parton distributions, 
which also predicts a sign difference for the quark Sivers function between the SIDIS and 
Drell-Yan processes. More complicated results have been found for the SSAs in the hadronic 
dijet-correlation 0, [s^, where a normal TMD factorization breaks down [31]. The reason is 
that the eikonal propagators from the initial and final state interactions in dijet-correlation 
process do contribute poles in the cross section [30|, l3l| • Because of this, the Ward identity 
is not applicable, and the standard TMD factorization breaks down, although a modified 
factorization may be valid if we modify the definition of the TM D p arton distributions to 
take into account all the initial and final state interaction effects [29]. In the fragmentation 
process, as we discussed in our paper, the eikonal propagators do not contribute to an 
imaginary part, and the Ward identity is applicable. We have shown this in our explicit 
calculations including one-gluon and two-gluon exchange contributions. 

There has been discussion about the twist-three quark-gluon correlation contribution in 
the fragmentation function, especially for the Collins effects 23|, |43|] . It will be interesting 
to further understand these contributions following the analysis in this paper, and discuss 
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35]. 



the universality issues in a more general ground 
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